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ABSTRACT 

Fermi observations of high-energy gamma-ray emission from GRB 0809 16C shows that its spectrum is 
consistent with the Band function from MeV to tens of GeV. Assuming one single emission mechanism dom- 
inates in the whole energy range, we show that this spectrum is consistent with synchrotron origin by shock- 
accelerated electrons. The simple electron inverse-Compton model and the hadronic model are found to be less 
viable. In the synchrotron scenario, the synchrotron self-Compton scattering is likely to be in the Klein-Nishina 
regime and therefore the resulting high-energy emission is subdominant, even though the magnetic field energy 
density is lower than that in relativistic electrons. The Klein-Nishina inverse-Compton cooling may also affect 
the low-energy electron number distribution and hence results in a low-energy synchrotron photon spectrum 
n(i^) oc zy"' below the peak energy. Under the framework of the electron synchrotron interpretation, we constrain 
the shock microphysical parameters and derive a lower limit of the upstream magnetic fields. The detection 
of synchrotron emission extending to about 70 GeV in the source frame in GRB 0809 16C favors the Bohm 
diffusive shock acceleration if the bulk Lorentz factor of the relativistic outflow is not significantly greater than 
thousands. 

Subject headings: gamma rays: bursts 



1. EVITRODUCTION 

It was recently reported that Fermi satellite has detected 
strong > lOOMeV emission from a very energetic long- 
duration burst GRB080916C (Abdo et al. 2009). At a red- 
shift of z = 4.35 ± 0.15 (Grenier et al. 2009), the burst is 
the most energetic one ever, with an isotropic gamma-ray en- 
ergy £^ ~ 8 X 10^"* ergs, which is released over a duration of 
Tgo ~ 60s. Equally remarkably, more than ten photons with 
energy above GeV are detected, with the highest energy one 
at 13 GeV (in the observer frame). The spectra of all five 
time intervals, designated as times a-e in the light curves of 
GRB080916C (Abdo et al. 2009), are wefl fit by the empir- 
ical Band function (Band et al. 1993), which smoothly joins 
low- and high-energy power laws. The high-energy power 
law extends to GeV energies, without any additional spec- 
tral component visible. The peak energy of the spectra dur- 
ing these intervals is around Sp ~ 400KeV- IMeV. Except 
during the first time interval, the low-energy and high-energy 
photon spectral indices of the prompt emission are constant 
and consistent with a ~ -1 .0 and /3 ~ -2.2 respectively. With 
its high temporal resolution, INTEGRAL detected the tem- 
poral variability of the KeV/MeV emission on time scales as 
short as 100 ms with high statistical significance (Greiner et 
al. 2009). So the variability timescale in the local source 
frame is r,, < 100ms/(l +z) = 20ms. 

The nonthermal synchrotron radiation by electrons has been 
suggested to be a possible mechanism for the 10 KeV-MeV 
emission (see Meszaros 2006 and Zhang 2007 for recent re- 
views), but one famous problem remains so far, i.e. the low- 
energy photon spectral index a is incompatible with the index 
-3/2 that is expected from fast-cooling electrons (e.g. Preece 
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et al. 1996; Ghisellini et al. 2000). Electron inverse Compton 
emission has been a competitive mechanism (e.g. Panaitescu 
& Meszaros 2000). The fact that one single spectral compo- 
nent fits the spectrum of the prompt emission from 10 KeV 
to GeV in GRB080916C suggests that one emission mecha- 
nism dominates in this whole energy range*. In this Letter, we 
study the constraint that this puts on the emission mechanism. 
Abdo et al. (2009) mentioned as one of the possibilities that 
the delay of high-energy gamma-ray emission relative to low- 
energy emission in GRB080916C could be a result of longer 
acceleration time needed for higher energy protons or nuclei 
in hadronic emission models. In accordance with this, we 
also study whether such hadronic models could be a possible 
mechanism that produces the KeV/MeV to GeV emission in 
GRB080916C. 

2. THE SYNCHROTRON MODEL AND PARAMETER CONSTRAINTS 

Assuming that in GRB shocks, fractions of and e^, of 
the shock internal energy are converted into the energy in the 
magnetic field and electrons, respectively. To ensure a high 
radiation efficiency for the prompt emission, it is usually as- 
sumed that the electrons are rapidly cooling, so the energy 
density in gamma-ray emission is equal to the electron en- 
ergy density Ug. The magnetic field is given by 
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where R is the radius of the shock, is the luminosity in 
gamma-ray emission and F is the bulk Lorentz factor The 
detection of GeV photons suggest that the emission region 
has a bulk Loentz factor F > 10^(Grenier et al. 2009; Abdo et 
al. 2009) ^. From the casuality constraint, the emission radius 
is R = T^cti,. In the synchrotron model for the 10 KeV-GeV 

^ Dropping the assumption of the same origin for both MeV and 
high energy emission, Li (2008) explained the high energy emission in 
GRB080916C as delayed, "residual" emission from subsequent collisions at 
larger and larger radii in the baryonic outflow. 

' The limit on the bulk Lorentz factor would be more stringent when we 
have an IC TeV component whose flux is above the synchrotron extension. 
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emission, by use of Sp = hvp = 'Pvj!^^llX, one can derive 
the Lorentz factor of electrons that radiate at the GRB peak 
energy Ep, 



_ .47mv;£ps 1/2/6^x1/4/ 2L^s-l/4 
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(2) 



where (f)iy ~ 0.5 is the coefficient defined in Wijers & Galama 
(1999) and q is the electron charge. Define 7t as the Lorentz 
factor of electrons below which the scattering with peak en- 
ergy photons are in the Thomson scattering regime, i.e. 

7r = rm,cVep = 250r3(£p/2MeVri. (3) 

Unless Eg < lO~'^eB, which is unreasonable in terms of the 
burst energetics, the IC scattering between 7„ electrons and 
the bulk of the gamma-ray emission should be in the Klein- 
Nishina (KN) regime. 

The KN Compton cooling of electrons may affect the low- 
energy electron distribution at 77 < 7 < 7,,, and hence affect 
the low-energy spectral slope below Sp (e.g. Rees 1967; Deri- 
shev et al. 2003), as we show below. Consider a population of 
electrons injected into a box with magnetic field B in a power 
law form N(j) oc 7"^ for 7 > 7m. These electrons will cool 
down rapidly through synchrotron and/or IC radiation^. As 
the electron cools, its energy changes as 7, so we have 

i^F,[l+kij)U^/UB]^J (4) 

for 7 < 7m, where v is the synchrotron frequency of 7- 
electrons and ^(7) accounts for the reduction of the effec- 
tive photon energy density for IC scattering of the 7-electrons 
due to the KN effect. Define Iii^kn = ^meC^f-y as the criti- 
cal energy of the photons with which the IC scattering of 7- 
electrons is just in the KN regime. For a low-energy photon 
spectrum of the form i/F/y v'^ ( v < Vp),-we. have 



(5) 



for 7r < 7 < 7m and k~\ for 7 < 77. In the case of 77 < 
7 < 7„, if W^/Ub ^ 1, i.e. the IC cooUng is still dominant 
even though the scatterings are in the KN regime, one can 
obtain vF^, oc 7*"'"' oc iy(''+iV2^ where we have used oc 7^ in 
the last step. Equating this synchrotron spectral index with 
the initially assumed spectral index Fv oc v°~^, one can derive 



5 = l,F„ocz/" 



(6) 



This can explain the observed low-energy photon index of a = 
-1 .0 in GRB080916C if the condition. Us < k(j)Uj for 7t < 
7 < 7m, is satisfied. 

On the other hand, the high-energy spectrum of 
GRB080916C above Sp is consistent with the synchrotron 
spectrum produced by fast-cooling electrons above 7„„ i.e. 
Fi, oc i/"''/^ with p = -2(1 + P2) = 2.4. The dominance of syn- 
chrotron cooling above Sp implies that Ub > k( j)U^ for elec- 
trons with Lorentz factor "f jm- Therefore, we find that, 
at 7 = 7m, Ub — kU-y. Since Uj ~ Ue for fast-cooling elec- 
trons, the requirement Ub — k{'y„)U-y translates into Ce/cB — 
ep/TmeC^, which gives 



7- 
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2MeV 



(7) 



as is the case in the synchrotron scenario. Using the formula in Lithwick & 
Sari (2001), however, one can find that the limit is only increased by a factor 
smaller than two. 

^ The synchrotron cooling time of 7m electrons is shorter than the dynamic 
time t' =R/rc = Tt^ as long as ee/es < 10*r3'*^\-i2(£p/2MeV)^/3L^,54. 



Note that the transition region between the two asymptotic 
power-laws at low and high energy ends in the Band func- 
tion is rather wide, so the above requirement, Ub — ki-fmWj, 
should be regarded as an order of magnitude of estimate. In 
addition, this requirement appUes only to large Sp bursts, be- 
cause for low Ep bursts, the IC scattering may be no longer 
in the Kelin-Nishina regime. A signature that high Ep bursts 
have a preferentially close to -1 can be seen in the analysis of 
Preece et al. (1996). For some low Ep bursts that have a ~ -1, 
some other mechanisms may be at work. 

For electrons with 7 < 7t, the IC scatterings with peak 
energy photons are in the Thomson scattering regime, so 
kij) = 1 and A^(7) oc 7"^, leading to a conventional fast- 
cooling photon spectrum of F^, oc Observations show a 
single power law spectrum oc from 10 keV to ^ MeV in 
GRB080916C, implying that 7m/7T > 10, and one can there- 
fore obtain a constraint 



^ > 1 »-2 / P \ 
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Due to that the IC scatterings between 7^ electrons and the 
peak energy photons with energy Ep are in the KN regime, the 
IC emission peaks at 



- ir^V/4r-l/4p2J/2/_eLxl/2j Y 
~ ^^ej S,54-^ 3V24MeV-' ^ 
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with a flux vFlf iE^i = hvf ) = i'(7m)z/pF,y„ ~ i^p/v„, where Y is 
the Compton parameter. For a flat synchrotron spectrum with 
/? ~ -2 above i^p, it is natural that the IC component is not 
seen at high energies since vF^{Eh = 70GeV) > vFl^iEh = 
70GeV) for GRB080916C. 

In the above, we have not assumed any model for the dissi- 
pation mechanism of the shocks. In the popular internal shock 
scenario, the typical Lorentz factors of the shocked electrons 
is 7m = ee{mp / me)T sh, where Tsh is the shock Lorentz factor, 
which is equal to the relative Lorentz factor of the two collid- 
ing sheUs. For GRB080916C, we have obtained a constraint 

7m ^ 5 X IQ^T^r^^^ltl^^i:^^). So if internal shock appUes 
to GRB080916C, we would need a large relative Lorentz fac- 
tor, Tsh - ^OeeT'T^L^'l^ltl'^i^). This could be caused by 
the interaction among the shells whose Lorentz factors have 
a large contrast (Yu et al. 2009). Of course, the shock could 
also arise from the magnetic reconnection or turbulence (e.g. 
Thompson 1994; Meszaros & Rees 1997; Lyutikov & Bland- 
ford 2003; Narayan & Kumar 2008; Zhang & Pe'er 2009) and 
we do not have the estimate of the shock Lorentz factor from 
the first principle. 

3. ALTERNATIVE MODELS FOR MEV- 10 GEV EMISSION? 

3.1. The one-zone SSC scenario 

Let's explore whether the simple IC scenario (i.e. one 
zone synchrotron self-Compton (SSC) scenario) can explain 
the single power- law spectrum from MeV to 10 GeV in 
GRB080916C. Suppose that the first-order SSC of electrons 
with energy 7,,, produce the peak emission Ep = 2MeV. Since 
this IC emission is not hidden by the synchrotron emission, 
one will expect that the 2nd order IC emission appears at 
high-energy if the 2nd-order IC peak is located within the 
observation energy window and that the 2nd-order IC scat- 
tering is still in the Thomson scattering regime. The fact 
that we did not see the 2nd-order IC component implies 
that 7^£p > 70GeV or 7m£p > TnieC^, so we have 7„ > 
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190(£p/2MeV)-'/2 or 7„ > 7t = 250r3(ep/2MeV)-i. Since 
Ep = hi'syn,p'ym, one obtains the synchrotron peak frequency at 
hvsyn,p = 55(7m/190)~^(£p/2MeV)eV. Then one can obtain an 
upper hmit of the magnetic filed B = hi'syn,p/{(t>v-^^^lm^) = 

140(7„/190)-^r3i(£p/2MeV)G. With this magnetic field, 
we can derive an upper limit for ee/e^.. In the case 
that 2nd-order IC is still in the Thomson regime, Ue = 
(F+ \)U^, so from es/e^ = UB/Ue ^ UbHYU^), we obtain 
F(eB/e«) = 3 x \0-\-fJ\90)-\ep/2MtWfL-%T\tl_^. By 
use of y = {te/esfl^ (Kobayashi et al. 2007), one get Y = 
\m^J\9Q)\ep/2MeSf)-'L'^%Tft-^\.^. As 7„ > 190, so the 
radiation energy in the 2nd-order IC will be E2nd,ic = YE^ > 
1.3 X lO^^erg, which is too large to be realistic. Such an 
energy crisis problem has also been found in the case of 
GRB080319B for the IC scenario of the prompt MeV emis- 
sion (Piran et al. 2008). 

On the other hand, if the 2nd-order IC is already in 
the deep KN regime (for 7m ^ = 250r3), Cs/ee = 
Us/Ue = Us/U-y. The 2nd-order Compton Y parameter 
is Y2ni = i<^e/<^By^^{lmhi:)~^ for a spectrum vF^ oc be- 
low Ep (see Eq.5). To get Find ^ 1> one needs 7^ > 
(ee/eB)'/'7T = 250(e,/eB)'/^r3(ep/2MeV)-i for ^ = 1. So 
from hvsyn,p = 55(7m/190)~^(£/,/2MeV)eV, one can obtain 
B < 45(ee/eBr^r3'(£p/2MeV)G. Combing this upper limit 
with the equipartition assumption in Eq.(l), one can get 

eel^B < Omv'^J^tll^^iep/lMtVfl^L^]!^], which is in conflict 

with the precondition ^ehB = Y^>l. This means that signifi- 
cant suppression of the 2nd-order IC emission by KN scatter- 
ing can not be fulfilled. So we conclude that the the simple 
one-zone SSC model does not work for the MeV to 10 GeV 
emission in GRB080916C. 

3.2. The hadronic scenario 

We first study whether the proton synchrotron emission 
can produce the MeV- 10 GeV emission of GRB080916C. 
The photon spectrum index above Sp, (3 = -2.2, implies that 
the proton distribution index is /? ~ 2.4 for fast-cooling pro- 
tons or a very steep index p ~ 3.4 for slow-cooling ones. 
In the proton synchrotron scenario, the observed peak emis- 
sion at Ep is produced by protons with a Lorentz factor of 

jp = (j^^^y^^- The synchrotron cooling time of these pro- 
tons in the comoving frame is t'^y^ = 6TrmpC/iaTm^jpB^). De- 
fine that the magnetic field energy density is a fraction of (_b 
of the comoving frame photon energy density, i.e. Ub = ^bU^. 
So the synchrotron radiation efficiency of the 7^ protons is 
T] = mm[l,t'j/t'^„], where = /?/rc is the dynamic time in the 
comoving frame, which is equal to the comoving frame vari- 
ability time, tj = t[, = Ff,,. As long as tl/t'^„<\, we have a 
radiation efficiency for 7^ protons 

= I; = 3 X io-*d/X?4Cf ri^(^)'/^ (10) 

Such a low radiation efficiency impUes an unreaUstically large 

energy in protons, a factor of l/?7 ~ 3 x 10^^^^''^ higher than 
the energy in gamma-rays. The radiation efficiency is quite 
low (?7 ~ 2 X 10"^) even for protons that produce the high- 
energy gamma-rays of energy > lOOMeV. So there is no 
room for the proton synchrotron model even in the assump- 
tion that the high-energy gamma-ray emission belongs to a 



different component than the MeV component. If the spec- 
trum above Sp is interpreted as arising from fast-cooling pro- 
tons, as required in the case of 2 < /? < 3, one would need 
£.B = Ub/U^ > 5 X 104L-ig4f2/3rf /'(2^)-2/3, which is also 
unreasonable, as the total energy in the magnetic field is too 
large to be reahstic for a GRB. 

Let's also explore the scenario of the secondary emission 
from hadronic photopion process. Detection of high-energy 
gamma-rays of energy greater than 10 GeV puts a constraint 
on the opacity of 77 absorption. As the hadronic pj opac- 
ity is related with the 77 opacity, a higher maximum photon 
energy, hence a lower 77 absorption opacity, would imply a 
lower hadronic radiation efficiency (e.g. Dermer et al. 2008). 
It is useful to express the hadronic 777 efficiency as a function 
of the the pair production optical depth t^-.,. Following Wax- 
man & Bahcall (1997), the optical depth for pair production of 

a photon of energy Sh is t^^(£/,) = = ri6 n«ll'h ' where 
/-.j -y is the mean free path. For the simplicity of calculation, we 
have assumed a photon spectrum P2 = -2 above Sp, which is a 
good approximation for GRB080916C. The fraction of energy 
lost by protons to pions is — YJ^^pj^peak for protons with 

energy greater than Bp = 6 x 10^^F^(£p/2MeV)"'eV (Wax- 
man & Bahcall 1997), where 10 ^**cm^ is the cross 
section of the pj reaction at the A resonance and ^peak — 0.2 
is the fraction of proton energy loss in one interaction. So the 
maximum photopion efficiency is 

/tt = 2 X 10 ^^li lOQi-y) '(iMeV-* '''t'I^^'i) (]]\ 

where Tj^ ~ (F/Fu^)"^ has been used in the last step and Fum 
is the minimum bulk Lorentz factor of the outflow required 
by the transparency of the high energy photon of energy Sh- 
As Tj-yish = 70GeV) < 1, we obtain an upper Umit of the pj 
efficiency, i.e. < 2 x 10"^ If the prompt MeV-10 GeV 
emission is interpreted as arising from the secondary emission 
of hadronic process, one would need an unreasonably large 
energy budget in relativistic protons. 

4. THE MAXIMUM SYNCHROTRON PHOTON ENERGY AND ITS 
IMPLICATIONS 

First we calculate the maximum synchrotron photon en- 
ergy that can be reached in GRB shocks. We assume that 
the relative motion between the upstream and downstream 
plasma is only mildly relativistic (such as in internal shocks). 
For an electron being shock accelerated, the residence time 
in downstream and upstream regions are respectively, tj = 
Kds'^ j/qBdC and f,', = nus'^^/qBuC, where e'^ is flie energy of 
accelerated electrons, Bj and B„ are respectively the magnetic 
fields in the shock downstream and upstream, and Kd,u ^ 1 pa- 
rameterizes the efficiency of shock acceleration, with Kd,u — 1 
corresponding to the fastest shock acceleration- the Bohm 
diffusive shock acceleration with the scattering mean free 
path equal to the particle gyroradius. It is generally assumed 
that the downstream magnetic field is close to the equipar- 
tition with the shock internal energy, while the value of up- 
stream magnetic field is less clear. As fi„ < B^, the total 
acceleration time is dominated by upstream residence time, 
so f/,^.^. ^ Ks'JqBuC. The maximum energy of accelerated 
electrons in each region is determined by equating the res- 
idence time with the shorter one of the cooling time and 
the available dynamic time, i.e. f^^^ = min{f^^^;,f^^,^}. The 
cooling time in the downstream and upstream are, respec- 
tively, = 3OTec/(4crr({/ij^„-FA;(7){/^)), where Vs^ and 
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Ub„ represent the magnetic field energy density in downstream 
and upstream respectively, and is the photon energy den- 
sity. In downstream region, Ub^ 3> kijMjWy, so the maxi- 
mum electron energy is jM.d = ( KjarBj ^^^^' upstream re- 
gion, the magnetic field energy density could be lower than 
k{-fM,u)Uj, and in this case, the maximum electron energy is 
)i/2, where k(-fM„) = (n/jM,.)''^- So 



_ / 3qB„ 



3qB„ 



-n)^/^. As Bd > B,i, the electrons radiate 



lM.li - (- „ „. 

more efficiently in the downstream and therefore the rele- 
vant maximum Lorentz factor with the observed radiation is 
7m = min[7M.,,,7M,rf]- Depending on which of 7m.,, and jmj 
is larger, we divide the discussion into two cases: 

i)The "fM.d S iM.ii case. The maximum synchrotron photon 
energy is 



synM 



1 



= 0.2294|^7^,r = 55( — jrjGeV, (12) 



which is only dependent of the bulk Lorentz factor T of the 
relativistic outflow (0.2294 is the coefficient quoted from Wi- 
jers & Galama 1999). If the observed highest energy photon 
is produced by synchrotron radiation, from hi^syn.M ^ £h, we 



obtain 



< 



o.sr^ ( ) \ (13) 

V70GeV/ ' ^ ' 

which favors the Bohm diffusive acceleration if the bulk 
Lorentz factor r3 < a few. 

From the precondition, jM.d ^ Jm.u, we obtain a lower limit 
of the upstream magnetic field in this case. 
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:500(|^)^/«..„«-^/V^(3^,r;'^/\e,/2MeV)-'/2f,7!//G 



(14) 



ii)For the case ^m.u ^ iM.d, the maximum synchrotron pho- 
ton energy is hi^synMlM.,,) = 0.2294||^7^„r. So from 
hvsyn.M{lM.i<) ^ £h, we obtain a lower limit of the upstream 
magnetic field. 
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0.2294x3^Srfr 
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3/4 - 1/2 
35 

,3/4 5/8 y.-4. 



■^7,54-^ 3 2MeV 



(15) 



Combining this lower limit with the precondition 7m.„ ^ 
iM.d, we find that Eq.(15) is applicable only when e;, < 
SSk^'TsGcV. In both cases, the GRB shells that produce the 



' and Eg ^ Bff, ^ r . For an initial magnetic filed 
lO'^G within a volume of radius of 10*- lO^cm, the 



prompt emission must have a pre-shock magnetic field greater 
than - 500G at a radius ofR ~ 3 x lO'^cmr^f,,,-!. If the field 
lines in the expanding shell are frozen and the width of the 
shell is constant, the components then vary with distance as 

Br oc r 
ofBo- 

above limit is larger than the Br component, but still within the 
Be or B^ component. Of course, the above limit is also con- 
sistent with the hypothesis that the upstream magnetic field 
is significantly amplified by the particle streaming instability 
(Bell 2004)^. Interestingly, the shock compressed upstream 
magnetic field, B ^ AVshBu ^ 1500G, is similar to the as- 
sumed equipartition magnetic field in downstream (i.e. Eq. 1), 
which means that the field compression due to the shock is 
enough to explain the downstream magnetic field. 

5. SUMMARY AND DISCUSSIONS 

The single-component spectrum of GRB080916C from 
MeV to GeV puts useful constraints on the emission mech- 
anism. We found that the synchrotron mechanism from rel- 
ativistic electrons is consistent with the observed spectrum, 
while the simple one-zone electron IC and hadronic models 
are less viable. In the synchrotron interpretation, the SSC 
emission is found to be in the KN scattering regime and as 
a consequence, the IC component is not visible at high ener- 
gies even though the magnetic field energy density is smaller 
than that in the relativist electrons, i.e. eg < e^, as obtained in 
our case. We also suggest a scenario in which such a KN IC 
emission dominated regime can explain the low energy pho- 
ton spectral index of GRB 0809 16C. 

The delay of high-energy gamma-ray emission relative to 
the low-energy emission in GRB080916C is still a mystery in 
the electron synchrotron scenario. It could be due to that the 
energy distribution slope p of electrons during the first time 
interval (time a) is rather steep so that the high-energy emis- 
sion is suppressed or that the emission region has not become 
transparent for high-energy gamma-rays at early times. 
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' Recently, Li & Waxman (2006) constrained the pre-shock magnetic fields 
of GRB afterglow shocks by synchrotron X-ray afterglows, which also im- 



plies that the pre-shock magnetic fields may be amplified. 
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